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SPECIFICATION 

VARIANCE REDUCTION TECHNIQUE FOR LARGE SCALE RISK MANAGEMENT 

Field of the Invention 

This invention relates to methods for accelerating the convergence of Monte Carlo 
simulation for pricing portfolios of loans, mortgages or related financial instruments. 

Description of Related Art 

Financial institutions may carry large portfoUos of mortgages and loans, numbering even 
in the miUions of financial instruments. To evaluate the value of such portfolios, a number of 
Monte Carlo simulations are commercially available such as in the Oracle OFSA/Risk Manager 
software package. This particular software package is currently used by many large banks for 
evaluating their portfolios of mortgages, loans and deposits in what is called Asset and Liability 
Management (ALM), Speed is a key concern for ALM packages for a number of reasons. 

For example, estimating the value of such a portfolio using a Monte Carlo simulation 
with a single representative instrument is very inaccurate; therefore all instruments (which may 
number in the millions) are often processed individually by the simulation before accumulating 
the results. In addition, prepayment assumptions can be specified in a very sophisticated manner 
(at least in the OFSA/Risk Manager), which slows down the computations. Because of these 
factors, prior art Monte Carlo simulations tended to minimize the number of scenarios used 
while maintaining a satisfactory degree of accuracy. 
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A particularly difficult financial problem involves the evaluation of a portfolio of fixed 
rate mortgages, which are perhaps the most comphcated of all banking instruments because of 
the path-dependency of amortization and prepayment. A simple Monte Carlo simulation of such 
a portfoUo may require many thousands of samples before the simulation converges on a result 
with a suitable degree of accuracy. A number of variance reduction techniques may be used to 
reduce the required number of samples/scenarios used in a Monte Carlo simulation, including 
numerical integration using low discrepancy sequences (LDS) or control variates. The latter 
technique has not proven as popular as LDS or antithetic variables for several reasons. For 
instance, these techniques require an estimation phase to extract information about the system, 
adding overhead to the calculations. If only a few mortgages are in the portfoho, this extra 
overhead can be comparable to the time required to process the system with a simple Monte 
Carlo simulation. 

These considerations do not apply, however, in evaluating a large portfolio. The 
computational overhead of a control variate technique is dwarfed by the sheer time of calculating 
the cash flows of a large portfoho. Also mortgages do not change very much with time. This 
means practically that one could estimate todays portfolio based on information from yesterday's 
portfolio, easing the estimation task of a control variate technique. Thus, there is a need in the art 
for improved control variate techniques that reduce the amount of scenarios required for a Monte 
Carlo simulation to converge to an answer or desired result. 

Summary of the Invention 

The invention provides in one aspect a control variate technique for Monte Carlo 
simulation of the value of a portfolio of financial instruments. The control variate is formed by a 
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linear combination of an estimated value for the portfolio, an estimated aggregate value the 
portfolio^ and an expected value of the estimated aggregate value. To derive the approximate 
simulated portfolio, a number of Monte Carlo samples are used to interpolate functions 
representing the approximate principal and approximate cashflow for the portfoUo. 

Other and further advantages of the invention will be further understood and appreciated 
by those skilled in the art by reference to the following specification, claims and drawings. 

Brief Description of the Drawings 

The accompanying drawings, which are incorporated in and constitute a part of the 
specification, illustrate preferred embodiments of the invention, and together with the 
description, serve to explain the principles of the invention. Similar reference characters denote 
corresponding features consistently throughout the drawings. 

Figure 1 is a flowchart of a simple Monte Carlo simulation. 

Figure 2 is a flowchart of a control variate technique according to one embodiment of the 
invention. 

Figure 3 is a graph comparing data obtained from a Monte Carlo simulation using a 
control variate technique according to one embodiment of the invention with data obtained form 
prior art Monte Carlo simulations. 

Fig. 4 is a diagram of a computer hardware system with which the present invention can 
be implemented. 

Fig. 5 is an additional diagram of computer hardware system with which the present 
invention can be implemented 
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Detailed Description 

The control variate technique of the invention may be used to calculate the market value, 
or price (at time zero) of a portfolio of mortgages. By definition the market value of the 

portfolio is the sum of the market values F/ of each mortgage. As will be shown herein, this 
market value may be represented by an expected value of a stochastic integral. The integrand 
will depend upon the (risk-fi-ee) interest rate as well as the financial characteristics of the 
particular mortgage such as its coupon rate and the starting principal. To estimate a value for the 
integral, the integrand is sampled over the payment times of the mortgage. Various stochastic 
models or term structures may be used to generate the interest rates at these payment times. In a 
simple Monte Carlo simulation, values or samples of the integrand are calculated using an 
interest rate scenario produced by sampHng a stochastic interest rate model with a pseudo- 
random number produced by a pseudo-random number generator. Multiple scenarios are 
produced and averaged to estimate the market value of the portfoho. 

A popular stochastic model or term structure for the risk-fi-ee interest rate r is the Hull 
and White one- factor model. The stochastic differential equation (SDE) for the Hull and White 
one-factor model is given by: 



where a, a are constants, ^ is one-dimensional Brownian motion in the risk-neutral measure, 
and <9 is a deterministic fixnction of time that optimally fits the current yield curve. This means 
that given K observable discount bond prices for A: = l,...,ir , a perfect choice of 0 

would result in: 



dr^ = (0^ - ar^)dt adW^ 



(1) 



P(T,)-E[e 




(2) 
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where E denotes time-zero expectation in the risk-neutral measure. By standard no-arbitrage 
arguments, the price of each mortgage then satisfies: 



where T-^ diXQ payment times, XhQ cash flows 0{T.^)^xt F^, measurable random variables that 

depend on mortgage characteristics^ and n is the number of payments times required to retire the 
mortgage. Consider the case of fixed rate mortgages with prepayment optionalities. For 
convenience of notation but without loss of generahty, it may be assumed that the maturity of 
each mortgage is 7 = 30 years, and that the initial principal of each mortgage of the pool is the 
same. However, a typical mortgage portfolio consists of mortgages issued at different dates in 
the past, resulting in distinct values for the outstanding principal at time zero of each 
mortgage. Moreover, the coupon rate is different for each mortgage. 

A Monte Carlo simulation is currently the most popular technique to calculate a value for 
equation (3). In that methodology, practitioners need to empirically estimate prepayment 
behavior. The most general endogenous model of prepayment is that of a fimctional of the rate, 

i.e. the time-/ prepayment rate is a functional : C([0,/],31) (8) 91^ ^ 9?^, where C{[0,t\SR) is 
the space of 9? -valued continuous functions of t . This prepayment model captures path- 
dependent behavior such as bumout and index lags. An example prepayment model is described 

further herein. The cash flow rate per unit of time and principal p follow then as: 



Vi = E[f^e^^ '''' a {T^] 



(3) 




+ r,{r,p\7t')pi ifpi>Q 



otherwise 



(4) 
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dpi 



dt 




(5) 



with pI = 7C^ . 

In a Monte Carlo simulation, the integral of equation (3) is sampled at discrete points to 
estimate its value. Thus, equations (1) and (4) must be discretized also. The Euler scheme may 
be used to discretize equation (1), Calculating a discrete cash flow from equation (4) is 
described by F. Fabozzi and F. Modigliani, Mortgage and Mortgage-backed Securities Markets, 
Harvard Business School Press, 1992. The total cash flow rate C is the sum of each mortgage 
cash flow rate and the total principal p is the sum of each mortgage principal. 

A flowchart for a simple, prior art Monte Carlo simulation of equation (3) is illustrated in 
Figure 1. After starting the simulation, a random number generator generates a pseudo-random 
number to initiate the calculation of an interest rate scenario at step 105. For example, if the 
simulation is being implemented in C++, the function random_fraction() will provide a random 
number in the closed interval [0, . . . 1 ] . Because a computer cannot generate a truly random 
number, the number is merely pseudo-random. With the appropriate weighting to maintain the 
required probabiUty density, the pseudo-random number is plugged into the stochastic interest 
rate model at to generate interest rate samples at the payment times of the mortgage being 
evaluated. The stochastic interest rate model or term structure may be a discretized version of 
equation (3). Altematively, other interest rate models may be used such as the Rendleman and 
Bartter model, the Vasicek model, or Hull and White two-factor model. Rather than use a 
pseudo-random number, LDS values could also be used in the stochastic interest rate model in 
step 105 to generate the interest rate samples such as described in U.S. Pat. No. 5,940,810, the 
contents of which are hereby incorporated by reference as if fully set forth herein. At step 110, 
the interest rate samples and financial characteristics of the mortgage, such as its starting 
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principal and coupon rate, are used to generate cash flow values at the payment times. These 
values are used to calculate samples of the integrand in equation (3) to estimate, for one scenario 
in the Monte Carlo simulation, the value for a particular mortgage at step 115. Steps 1 10 and 
1 15 are repeated at step 120 until a value has been estimated for all the financial instruments in 
the portfolio. 

The sum of all the values of all the mortgages calculated in step 120 provides an estimate 
of the value of the portfolio for one particular scenario in the Monte Carlo simulation. Many 
such scenarios must be repeated at step 125 depending upon the variance of the results. In a 
large portfolio, many different scenarios must be calculated before the simulation converges on a 
suitably accurate result. Because the process illustrated in Figure 1 implements no variance 
reduction techniques, it may be denoted as a simple or crude Monte Carlo simulation. 

A control variate technique reduces the number of scenarios required for a Monte Carlo 
simulation to converge upon a result. The Monte Carlo simulation discussed with respect to 
Figure 1 may be generalized as providing an estimate of a fimction E(g(X)), where X is a random 
vector representing the interest rate samples at the payment times (a single Monte Carlo 
scenario). If another function f is derived whose expected value of f(X) is known, the 
convergence of the Monte Carlo simulation may be accelerated (assuming appropriate selection 
of a suitable f(X)) by calculating samples of a new function W: 

W-(g(X)-f(X) + m) (6) 
where m is the expected value of f(X). The function W may be denoted as the control variate. 
The inventor has discovered a control variate that converges to a result more than two orders of 
magnitude faster than a simple Monte Carlo simulation. 
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In the control variate technique of the invention, the function f(X) represents the value of 
an aggregated approximate portfoho. Note that it is much less time-consuming to simulate the 
approximate portfoho than the true portfolio because the former is aggregated. As used herein, 
when references are made to an "aggregated" value such as an "aggregated cash flow value," it 
will be understood that this refers to a value of the approximate portfoho, which does not 
simulate or model the financial instruments in the portfoho individually. To simulate the 

approximate portfolio, a functional form for its cash flow C- and principal for each month 
/ (payment time) may be interpolated. The inventor has discovered that using two state 
variables, the interest rate and the approximate principal itself, provides a powerful control 
variate technique. The model of the approximate portfolio may be given by: 

q^lLi{r{iM\p,,m) (7) 

p,^,=:v{r{iAtlp.JAt) (8) 
The functions ju and v are determined by the information extracted in a first set of Monte Carlo 
scenarios, e.g., co = 1,...,50, using steps 105 through 125 in Figure 1. Because these scenarios 
are used to model an aggregate portfolio, they may be denoted as "adaptive" scenarios. Note that 
the value of 50 scenarios is arbitrary. The number of scenarios required to approximate the 
functions ju and o is substantially less, however, than the many thousands of scenarios required 
in a simple Monte Carlo simulation. 

From equation (3), the control variate W is thus given by: 

Z-l rMS . * — 1 -I ras ^ ^ — » - r.as ^ , 

e'' C^(7;.0-Ze^» C,(«)-J;£[e CJ (9) 

where C-(6)) is given by equation (7) and represents the cash flow for the approximate portfolio 
in the ith month (or ith payment time). An examination of equation (7) shows the dependence of 
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C. {(X)) on the two state variables in the function // : the interest rate for the ith month, r(zA^, (o) , 



and the principal for the approximate portfoUo at the ith month, p^^) . Interpolating this 
function directly from the two state variables requires an interpolation in 3 dimensions. While 
such an interpolation can be done, there is a danger of overfitting the data in such a case. 

Rather than interpolating in three dimensions, an interpolation in just two dimensions can 
be accomphshed in the following manner. To decouple the dependence on the two state 
variables, the aggregate portfoho's cashflow for the ith month from only the payments (denoted 

as Cf""^ ) is distinguished from the cashflow resulting from prepayments, and the dependence of 
Cf^^ upon the principal is interpolated. In addition, the average prepayment rate for the 
aggregate portfoho in the ith month (denoted as f- ) is derived as a function of r(zAO . Finally, 
the principal for the {ith + 1) month resulting only from the average prepayments in the ith 
month (denoted as ^ ) is distinguished from the change in principal resulting from payments in 

the ith month, and the dependence of upon the principal is interpolated. 

The adaptive scenarios are run (steps 105 through 125 as shown in Figure 1) and the 
following statistics are extracted, for each time il^t and scenario co : 



where 1 is the indicator function, and j is the index for the jth financial instrument. A functional 
form for as a function of r{iAt,6)) may be derived in the following fashion. After graphing 
the average prepayment in ascending order of r(iAt, co) , a decreasing piecewise Hnear function 




(10) 




(11) 
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(with, e.g., 5 evenly spaced anchor points) is interpolated from the data to represent f .(r) , where 
f^{r) represents the prepayment rate for the aggregate portfoho. Similarly, the total payment is 
graphed in ascending order of Pi{o)) , to which an increasing piecewise linear function (with, 

e.g., five evenly spaced anchor points) C-''^{p^) is fitted. Finally, the ratio of p-^^ over (one 
minus the average prepayment for the ith month) is plotted in ascending order of p. to linearly 
interpolate the function , The function // and o are then defined by: 
Kr.pm = Cr{p)^pr,{r) (12) 

o{r,p4l^t)^pr{p){\''r,{r)) (13) 
It will be understood that the use of 5 anchor points is arbitrary as well as using piecewise linear 
functions. Other functional forms could be used with varying number of anchor points to 

approximate the functions f .(r) , Cf"^ , and p"!^^ . 

An additional set of Monte Carlo scenarios may be used to form the control variate W of 
equation (9). A flowchart for the control variate technique of the invention is illustrated in 
Figure 2. Upon starting the simulation, the statistics necessary to derive the average prepa5m[ient 
and total payment (equations (10) and (1 1)) are generated in a set of adaptive Monte Carlo 
scenarios at step 200. These scenarios are generated in an analogous fashion to steps 105 
through 125 described with respect to Figure 1 , The average prepayment and total prepayment 
may then be calculated using discrete forms of equations (10) and (1 1) at step 205. At step 210, 

the functions y^ir), Cj'''^ (p^) , mid pf^\p') are interpolated. 

To begin calculating the control variate W given by equation (9), the simulation begins 
another set of Monte Carlo scenarios at step 215 by generating an interest rate scenario at step 
215. At step 220, the interest rate scenario is used to calculate a value for the first term in 
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equation (9), This gives a value of the simulated portfolio according to one Monte Carlo 
simulation as described in steps 110, 1 15, and 120 with respect to Figure 1. The value of the 
simulated aggregated portfolio may now be calculated based upon the interest rate scenario 
provided by step 215. Note that equation (10) predicts the value for the aggregate principal for 
the {ith +1) month based upon the interest rate and the aggregate principal for the ith month. 
Thus, based upon the known starting principal values for the mortgages in the portfoUo and the 
interest rate samples, the aggregate principal values may be sequentially derived using equation 
(12) at step 225. In turn, the interest rate scenario and the aggregate principal values are 
substituted into equation (11) to calculate the aggregate cash flow values at step 230. These 
values are used in a discretized approximation of the second term in equation (9) to give a value 
of the simulated aggregate portfoUo according to one Monte Carlo simulation at step 235. Not 
illustrated in the flowchart is a step to calculate the expected value of the simulated aggregate 
portfolio. As will be described further herein, this term may be calculated by a simple Monte 
Carlo simulation or by a lattice approximation. The resulting expected value and the values from 
steps 235 and 220 are linearly combined at step 240 to calculate the control variate. Additional 
Monte Carlo scenarios are repeated at step 245 until the simulation converges to a result at step 
250. 

The aggregate cash flow and principal given by equations (7) and (8) is a rough 
approximation of the values given by equation (4)-(5). For instance, it does not capture burnout 
or index lags. The inventor has discovered, however, that most of the path-dependency of 
residential mortgages comes from principal amortization. Thus, an inventive solution is to take 
the approximate principal as a second state variable. This is a relatively optimum selection of 
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state variables since taking more than two state variables complicates the estimation phase, 
requiring interpolation in more than two dimensions, 

A simple Monte Carlo simulation may be used to calculate (with a very large number of 
scenarios) the expected value term for the control variate since the value of the aggregate 
portfolio is inexpensive to compute. However, computing the expected value with a lattice, i.e. a 
sparse representation of a Markov chain, provides a hard bound on the bias so that this technique 
is preferable. This is possible because the two state variables {r,p) constitute a joint Markov 
process. 

The lattice approximation of the invention is three dimensional and is an extension of the 
two-dimensional lattice discussed in Hull, Options, Futures, and Other Derivatives, p.438, 
Prentice-Hall, 1997, 3'^ Ed. to model instruments whose cash flows C depend on the interest 
rate only. The main features of the latter are now discussed to provide backgroimd for the novel 
lattice approximation of the invention. In Hull and White's two-dimensional lattice, time takes 
the discrete values ilst , and the discretized interest rate process r^{i) takes the discrete values 
a^ j = E[r{i/S.t)] + j/Sr , where the rate-step Ar is appropriately chosen for error minimization. 
The transition probabiUties at node {ij) of the lattice are zero, except for 3 
nodes (r + 1, - 1) , (z -f 1, y^^J , and {i + IJ^^^ + 1) ; they are such that the first two moments of 
the discretized process r^ for the rate equal the first two moments of the process r over a time- 
step . To calculate the market value V , equation (3) is discretized: 



where Q{iJ) is the price of an Arrow-Debreu security that pays 1 if node {iJ) is reached, i.e.: 



n 



(14) 



1=1 j 



Q(i,j) = E[e 



(15) 
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Note that this lattice has only two dimensions: time (through i) and interest rate (through y). In 
an embodiment of the lattice used to calculate the expected value term for the control variate, a 
dimension is added to account for dependence on the principal, such that the node {ij,k) is a 
time ibJ: state where the rate equals a^ j as before, and the principal equals k/Sp for some value 
of the principal step Ap . Conditional on that information, the principal at time {i + 1)A^ will be 
equal to o(a^ j,kApJAt) because of equation (8). At each node {i,j\k), k^^^ can be determined 
by rounding, so that: 

oia-j,kApJAt) -dp/2< k^^,dp < o(a-j,kApJAt) + dp/2 (16) 
Since the dynamics of the principal are deterministic between two time-steps, it is only possible 
to jump from node {ij\k) to nodes (/ + IJ^^, + hKext) ^ + 1' 

(} + ^ J next ~ 1' ^next) ' ^hc trausltion probabilities are the same as in the two-dimensional lattice. 
The market value of the approximate portfolio is then: 

Z Z Z 20'. J\ k)juU^. kAp, I At) ( 1 7) 

where QiiJ.k) isthepriceof an Arrow-Debreu security that pays 1 if the node {ij,k) is 
reached. 

Numerical results. An embodiment of the control variate technique of the invention was tested 

on a portfoUo of 100 mortgages. The test used pi - ; , p' = 0.062 + 0.0002 j - 0.02 * 

(77100)^for j = 1,...,100. The rate spread x and the prepayment function y are chosen to be: 

, ; r(iAt) + r((i~l)At) ^^^^ 
X = max[/?^ - -^^ — ' — J ^—^ " 0,005, 0] (18) 

fmin[0.01 + 10x + 5000x', I]/ At ifz<4; 
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The term structure parameters area = 0.1, cr = 0.014, and the yields of the discount bonds are 
given in the following Table 1. Cubic splines are used to interpolate these yields in order to fit 6 
according to equation (2). 



Table. 1. Original Yield Curve 


Maturity 


0.08 0.17 0.25 0.33 0.42 1 


Yield (in %) 


5 5.23 6 6.2 6.5 6.5 



The results are shown in Figure 3. In this test, the interest rate samples were generated using 
LDS samples for input into a stochastic interest rate model rather than with a pseudo-random 
number generator. Figure 3 compares the absolute value of the relative error (compared to a 
simulation of 100,000 pseudo-random numbers) of three estimators of the price of the mortgage 
portfolio: LDS, LDS with control variates, and pseudo-random numbers. For the latter, the 
average absolute value of the integration error is shown over 10 different seeds. For control 
variates, the 50 scenarios from the adaptive scenarios were not shovm in the results since in 
practice this data is readily available from previous simulations of the portfolio. As mentioned 
earlier, the overhead to generate the control variate is negligible. This experiment demonstrates 
that, for a mortgage portfolio, integration error is about 5 times lower with LDS than with 
pseudo-random numbers. If control variates with LDS is used, it is about 20 times lower. The 
variance of the pseudo-random estimator is between 100 and 300 times higher than the variance 
of the control variate estimator with LDS. 

An example computer system for implementing an embodiment of the invention 
is illustrated in Figure 4. Referring to Fig. 4, in an embodiment, a computer system 720 includes 
a host computer 722 connected to a plurality of individual user stations 724. In an embodiment, 
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the user stations 724 each comprise suitable data terminals, for example, but not Hmited to, e.g., 
personal computers, portable laptop computers, or personal data assistants ("PDAs"), which can 
store and independently run one or more appHcations, i.e., programs. For purposes of 
illustration, some of the user stations 724 are connected to the host computer 722 via a local area 
network ("LAN") 726. Other user stations 724 are remotely connected to the host computer 722 
via apubhc telephone switched network ("PSTN") 728 and/or a wireless network 730, 

In an embodiment, the host computer 722 operates in conjunction with a data storage 
system 731, wherein the data storage system 731 contains a database 732 for storing financial 
information about a portfolio of financial instruments that is readily accessible by the host 
computer 722. Note that a multiple tier architecture can be employed to connect user stations 
724 to a database 732, utilizing for example, a middle application tier (not shown). 

In altemative embodiments, the database 732 may be resident on the host computer, 
stored, e.g., in the host computer's ROM, PROM, EPROM, or any other memory chip, and/or its 
hard disk. In yet altemative embodiments, the database 732 may be read by the host computer 
722 fi-om one or more floppy disks, flexible disks, magnetic tapes, any other magnetic medium, 
CD-ROMs, any other optical medium, punchcards, papertape, or any other physical medium 
with patterns of holes, or any other medium fi'om which a computer can read. In an altemative 
embodiment, the host computer 722 can access two or more databases 732, stored in a variety of 
mediums, as previously discussed. 

Referring to Fig. 5, in an embodiment, each user station 724 and the host computer 722, 
each referred to generally as a processing unit, embodies a general architecture 805. A 
processing unit includes a bus 806 or other communication mechanism for communicating 
instructions, messages and data, collectively, information, and one or more processors 807 
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coupled with the bus 806 for processing information. A processing unit also includes a main 
memory 808, such as a random access memory (RAM) or other dynamic storage device, coupled 
to the bus 806 for storing dynamic data and instructions to be executed by the processor(s) 807. 
The main memory 808 also may be used for storing temporary data, i.e., variables, or other 
intermediate information during execution of instructions by the processor(s) 807. 

A processing unit may further include a read only memory (ROM) 809 or other static 
storage device coupled to the bus 806 for storing static data and instructions for the processor(s) 
807. A storage device 810, such as a magnetic disk or optical disk, may also be provided and 
coupled to the bus 806 for storing data and instructions for the processor(s) 807. 

A processing unit may be coupled via the bus 806 to a display device 811, such as, but 
not limited to, a cathode ray tube (CRT), for displaying information to a user. An input device 
812, including alphanumeric and other keys, is coupled to the bus 806 for communicating 
information and command selections to the processor(s) 807. Another type of user input device 
may include a cursor control 813, such as, but not limited to, a mouse, a trackball, a fingerpad, or 
cursor direction keys, for communicating direction information and conamand selections to the 
processor(s) 807 and for controlling cursor movement on the display 811. 

According to one embodiment of the invention, the individual processing units perform 
specific operations by their respective processor(s) 807 executing one or more sequences of one 
or more instructions contained in the main memory 808, wherein these sequences of instructions 
carry out the control variate technique of the invention. Such instructions may be read into the 
main memory 808 fi-om another computer-usable medium, such as the ROM 809 or the storage 
device 810. Execution of the sequences of instructions contained in the main memory 808 
causes the processor(s) 807 to perform the processes described herein. In alternative 
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embodiments, hard-wired circuitry may be used in place of or in combination with software 
instructions to implement the invention. Thus, embodiments of the invention are not limited to 
any specific combination of hardware circuitry and/or software. 

The term "computer-usable medium," as used herein, refers to any medium that provides 
information or is usable by the processor(s) 807. Such a medium may take many forms, 
including, but not limited to, non-volatile, volatile and transmission media. Non-volatile media, 
i.e., media that can retain information in the absence of power, includes the ROM 809. Volatile 
media, i.e., media that can not retain information in the absence of power, includes the main 
memory 808. Transmission media includes coaxial cables, copper wire and fiber optics, 
including the wires that comprise the bus 806. Transmission media can also take the form of 
carrier waves; i.e., electromagnetic waves that can be modulated, as in frequency, amplitude or 
phase, to transmit information signals. Additionally, transmission media can take the form of 
acoustic or Hght waves, such as those generated during radio wave and infrared data 
communications. 

Common forms of computer-usable media include, for example: a floppy disk, flexible 
disk, hard disk, magnetic tape, any other magnetic medium, CD-ROM, any other optical 
medium, punchcards, papertape, any other physical medium with pattems of holes, RAM, ROM, 
PROM (i.e., programmable read only memory), EPROM (i.e., erasable programmable read only 
memory), including FLASH-EPROM, any other memory chip or cartridge, carrier waves, or any 
other medium from which a processor 807 can retrieve information. 

Various forms of computer-usable media may be involved in providing one or more 
sequences of one or more instructions to the processor(s) 807 for execution. For example, the 
instructions may initially be provided on a magnetic disk of a remote computer (not shown). The 
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remote computer may load the instructions into its dynamic memory and then transit them over a 
telephone Une, using a modem. A modem local to the processing unit may receive the 
instructions on a telephone line and use an infrared transmitter to convert the instruction signals 
transmitted over the telephone Hne to corresponding infrared signals. An infrared detector (not 
shown) coupled to the bus 806 may receive the infrared signals and place the instructions therein 
on the bus 806. The bus 806 may carry the instructions to the main memory 808, from which the 
processor(s) 807 thereafter retrieves and executes the instructions. The instructions received by 
the main memory 808 may optionally be stored on the storage device 810, either before or after 
their execution by the processor(s) 807. 

Each processing unit may also include a communication interface 814 coupled to the bus 
806. The communication interface 814 provides two-way communication between the 
respective user stations 824 and the host computer 822. The conmiunication interface 814 of a 
respective processing unit transmits and receives electrical, electromagnetic or optical signals 
that include data streams representing various types of information, including instructions, 
messages and data. 

A communication link 815 links a respective user station 824 and a host computer 822. 
The communication link 815 may be a LAN 726, in which case the conmiunication interface 814 
may be a LAN card. Alternatively, the communication link 815 may be a PSTN 728, in which 
case the communication interface 814 may be an integrated services digital network (ISDN) card 
or a modem. Also, as a fiirther alternative, the communication link 815 may be a wireless 
network 730. 

A processing unit may transmit and receive messages, data, and instructions, including 
program, i.e., appUcation, code, through its respective conmiunication link 815 and 
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communication interface 814. Received program code may be executed by the respective 
processor(s) 807 as it is received, and/or stored in the storage device 810, or other associated 
non-volatile media, for later execution. In this manner, a processing unit may receive messages, 
data and/or program code in the form of a carrier wave. 

Those of ordinary skill will appreciate that many different embodiments of the invention 
may be constructed. Thus, it is understood that other embodiments of the present invention can 
be fabricated and be within the spirit and scope of the appended claims. 
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What is claimed is: 

1 . A method for financial estimation, comprising: 

(a) providing a portfolio of financial instruments having a schedule of payment times; 

(b) generating a plurality of interest rate scenarios by Monte Carlo simulation using a 
stochastic term structure; 

(c) calculating, for each fmancial instrument; a set of financial outcomes using a subset of 

interest rate scenarios; 

(d) interpolating, from the sets of fmancial outcomes, a first fimction that estimates the 
aggregate value of the portfolio when sampled at the payment times with an aggregate principal 
of the portfoho and a first interest rate scenario; 

(e) providing a second function that estimates a value for a fmancial instrument firom the 
portfolio when sampled at the payment times using a second rate scenario; 

(f) using an interest rate scenario fi-om the plurality of interest rate scenarios, calculating 
a value of the second function for each financial instrument in the portfolio and a value of the 
first function to estimate, respectively, a value for the portfoho and an aggregate value for the 
portfolio, and; 

(g) forming a control variate base upon the estunated value of the portfoho, the estimated 
aggregate value of the portfoho, and an expected value of the aggregate value of the portfoho. 

2. The method of claim 1 , further comprising: 

repeating steps (e) through (g) and averaging the resulting control variates to estimate a 
value for the portfoho. 
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3 . The method of claim 2, wherein the interest rate scenarios used in the repetition of steps 
(e) through (f) do not belong to the subset of interest rate scenarios used in step (c). 

4. The method of claim 3, wherein step (d) comprises deriving an average prepayment for 
the portfoUo as a fimction of the interest rate, deriving a third function giving the cashflow not 
resulting from prepayments for the portfoho as a function of the principal of the portfoUo, and 
deriving a fourth function giving the principal of the portfolio at an (i + l)th payment time based 
upon the change in principal predicted by the value of the average prepayment for the portfoUo at 
an ith payment time. 

5. The method of claim 4, wherein the aggregate principal at the (i+l)th payment time is 
given as the value of the fovirth function minus the product of fourth function times the average 
prepayment at the ith payment time. 

6. The method of claim 5, wherein an aggregate cashflow at the ith payment time is given as 
a sum of the third function and the product of aggregate principal times the average prepayment 
at the ith payment time. 

7. The method of claim 6, wherein the first fimction is a fimction of the aggregate cashflow. 
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8. A method for financial estimation of a portfolio of financial instruments, comprising: 

(a) providing a first function giving an aggregate value for the portfolio; 

(b) providing a second fimction giving a value for a financial instrument within the 
portfolio; 

(c) using an interest rate scenario, calculating a value for the portfolio using the second 
function and the aggregate value for the portfolio using the first function; 

(d) forming a control variate based upon the value for the portfoUo, the aggregate value 
for the portfolio, and an expected value of the aggregate value for the portfolio. 

9. The method of claim 8, wherein the first function gives an aggregate value for the 
portfolio as a function of an aggregate principal of the portfoho and an interest rate scenario. 

10. A computer program product that includes a computer-usable medium, the medium 
having stored thereon a sequence of instructions which, when executed by a processor, causes 
the processor to execute a process for financial estimation, said process comprising: 

(a) providing a portfolio of financial instruments having a schedule of payment times; 

(b) generating a pluraUty of interest rate scenarios by Monte Carlo simulation using a 
stochastic term structure; 

(c) calculating, for each financial instrument; a set of financial outcomes using a subset of 

interest rate scenarios; 

(d) interpolating, fi-om the sets of financial outcomes, a first function that estimates the 
aggregate value of the portfolio when sampled at the payment times with an aggregate principal 
of the portfolio and a first interest rate scenario; 
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(e) providing a second function that estimates a value for a financial instrument from the 
portfolio when sampled at the payment times using a second interest rate scenario; 

(f) using an interest rate scenario from the plurality of interest rate scenarios, calculating 
a value of the second function for each financial instrument in the portfoUo and a value of the 
first fiinction to estimate, respectively, a value for the portfoho and an aggregate value for the 
portfolio, and; 

(g) forming a control variate base upon the estimated value of the portfolio, the estimated 
aggregate value of the portfoUo, and an expected value of the aggregate value of the portfolio. 

1 1 . The method of claim 1 0, fiirther comprising: 

repeating steps (e) through (g) and averaging the resulting control variates to estimate a 
value for the portfoUo. 

12. The method of claim 1 1, wherein the interest rate scenarios used in the repetition of steps 
(e) through (f) do not belong to the subset of interest rate scenarios used in step (c). 

13. The method of claim 1 2, wherein step (d) comprises deriving an average prepayment for 
the portfoUo as a function of the interest rate, deriving a third fimction giving the cashflow not 
resulting from prepayments for the portfoUo as a fimction of the principal of the portfoUo, and 
deriving a fourth fimction giving the principal of the portfoUo at an (i + l)th payment time based 
upon the change in principal predicted by the value of the average prepayment for the portfolio at 
an ith payment time. 
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14. The method of claim 13, wherein the aggregate principal at the (i+l)th payment time is 
given as the value of the fourth function minus the product of fourth function times the average 
prepayment at the ith payment time. 

15. The method of claim 14, wherein an aggregate cashflow at the ith payment time is given 
as a sum of the third function and the product of aggregate principal times the average 
prepayment at the ith payment time. 

16. The method of claim 1 5 , wherein the first function is a function of the aggregate 
cashflow. 
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ABSTRACT 

A control variate Monte Carlo simulation for estimating the value of a portfolio of 
financial instruments is provided. The control variate uses a function dependent upon two state 
variables, the aggregate principal and the interest rate, to simulate the aggregate value of the 
portfolio. 
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